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CN) , We investigate the decay of the global entanglement, due to decoherence, of multiqubit systems 

interacting with a reservoir in a non Markovian regime. We assume that during the decoherence 
process each qubit of the system interacts with its own, independent environment. Most previous 
s ! \ works on this problem focused on particular initial states or families of initial states amenable of 

C*) ■ analytical treatment. Here we determine numerically the typical, average behaviour of the system 

^ ■ corresponding to random initial pure states uniformly distributed (in the whole Hilbert space of 

| n-qubit pure states) according to the Haar measure. We study systems consisting of 3, 4, 5, and 

6 qubits. In each case we consider also the entanglement dynamics corresponding to important 
particular initial states, such as the GHZ states or multiqubit states maximizing the global entan- 
glement, and determine in which cases any of these states is representative of the average behaviour 
associated with general initial states. 
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I. INTRODUCTION 

>; 

Decoherence is a quantum phenomenon that plays an important role in connection both with the foundations 
of quantum physics and with its technological applications |ln3jj. Decoherence is inextricably linked to another 
basic ingredient of the quantum world: quantum entanglement [J-0|. Indeed, the diverse effects associated with the 
phenomenon of decoherence are due to the development of entanglement between the system under consideration 
and its environment. This entanglement, arising from the interaction between an imperfectly isolated system and its 
surroundings, leads to the gradual disappearance of various quantum features exhibited by the system. These effects 
are at the core of the nowadays orthodox, decoherence-based explanation of the quantum-to-classical transition [l[ . 
• • ■ One of the most important consequences of decoherence is that the (internal) entanglement between the different 
parts of a composite system undergoing decoherence tends to decay. This manifestation of decoherence constitutes one 
■ of the main hurdles that has to be overcome for the implementation of quantum technologies that require entangled 
states as a resource @- The decay of entanglement under decoherence has been the focus of an intense research 
activity in recent years [8H20l]. A particularly interesting feature of this process, first pointed out by Zyczkowski and 
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the Horodeckis in a pioneering effort almost ten years ago [21], is that sometimes entanglement vanishes completely 
in a finite time. This remarkable phenomenon, known as entanglement sudden death (ESD), has been the subject of 
several theoretical and experimental studies [22i - [37| . 

Studies of ESD in Markovian scenarios suggested that, once entanglement has vanished due to decoherence, this 
valuable resource is irremediably lost. However, the decoherence process in non-Markovian regimes sometimes gives 
rise to an interesting new effect: entanglement sudden revival. The sudden revival of entan glem ent after a finite 
time interval with zero entanglement following an ESD event has been recently investigated in [271 28] . These works 
considered a system consisting of two independent qubits interacting with an environment at zero temperature in a 
non-Markovian regime. The entanglement dynamics of open multiqubit systems in non-Markovian regimes has also 
been considered in various recent works l29l. l30l, [3pl . In particular, multiqubit systems evolving from initial GHZ 
or W states have been studied in detail in [29l l30j|. On the other hand, the average, typical entanglement-related 
features corresponding to some families of initial states of two-qubit, three-qubit, and four-qubit systems have been 
investigated in 32]. 

The phenomenon of entanglement sudden revival allowed by non-Markovian regimes deserves detailed scrutiny: it 
is of clear interest from the fundamental point of view and may also have technological implications. The aim of the 
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present contribution is to explore (statistically) the typical behaviour of the entanglement exhibited by multipartite 
systems interacting with the above mentioned environment. In contrast to previous studies of entanglement evolution 
in non-Markovian regimes, that focused on particular states or particular families of states, we will explore the average, 
typical entanglement dynamics associated with random, general initial pure states uniformly distributed according 
to the Haar measure. We are also going to consider the average behaviour corresponding to mixed initial states 
consisting of a mixture of a (random) pure state and the totally mixed state. 

The article is organized as follows. A brief description of the system studied here is provided in Section II. Our 
main findings concerning the typical entanglement evolution of n-qubit systems are detailed in Section III. Finally, 
some conclusions are drawn in Section IV. 



II. SYSTEM-RESERVOIR MODEL 



We are going to explore the typical, average features characterizing the entanglement dynamics of an open system 
consisting of n independent qubits evolving in a non-Markovian regime. These n-qubits evolve from an initial en- 
tangled state generated through some previous, unspecified, interacting process. Each of the above mentioned qubits 
experiences the same decoherence process, which is due to the interaction of the qubit with its own environment. The 
evolution of the density matrix pi of the i-qubit is given by a completely positive trace-preserving map that admits a 
Kraus representation of the form, 



M 
3=1 

where Ej j — 1, . . . , M constitute an appropriate set of Kraus operators yielding a complete description of the single 
qubit's dynamics. The evolution of the total density matrix p corresponding to the entire n-qubit system can also be 
expressed in terms of the above Kraus operators, 



p(0) —>p(t) = J2 E n®---® E jn P(°) [En® ...®E jn ]l (2) 



We are going to assume that the dynamics of each qubit of the system is determined by a "qubit + reservoir'' 
interaction given by the paradigmatic Hamiltonian (27j . 



H = oj (t + (7- + uikblfik + ^.gfe(cr+6fe + £r_b£), (3) 
fe fe 

where loq is the transition frequency of an atomic two level system (qubit), <r+, <r_ are the atomic raising and lowering 
operators, bt, bk are the creation and annihilation operators of the reservoir's (electromagnetic) fc-mode, and u>k, 
<7fc denote the corresponding frequency and coupling constant (with the qubit). A possible physical realization of a 
system governed by the Hamiltonian ([3]) is given by a two-level atom interacting with a reservoir consisting of the 
electromagnetic modes of a high-Q cavity. Following j2?| we adopt for the reservoir an effective spectral density of 
the form 



ii \ 1 7o A 2 

JH = 2 ^-0,)' + A' " (4) 

The parameter A, defining the spectral width of the coupling, is then connected to the reservoir correlation time tj, 
trough Tfc ~ A -1 . On the other hand, the parameter 70 is associated with the system's relaxation time tr according 
to the relation tr ~ (see (2?| for details). Non-Markovian dynamics is obtained in the strong-coupling regime 
given by the inequality 70 > A/2. The time dependent density matrix describing the behaviour of a single qubit has 
the form 



P(t) = (poo(Q)+/»n(Q)(l-J>t))|0)(0| + poi(0)vS|0)<l| + p 10 (0)^|l)<0| + Pll (0)p t \l)(l\, 
where the (0) denote the elements of the initial density matrix, and pt is given by 



(5) 
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where d = y / 2'foX — A 2 [27| • In the present work we are going to focus on the non-Markovian regime, setting A = O.OI70 
in all our computations. 

It can be verified, after some algebra, that the evolution of the single qubit reduced density matrix elements given 
by ^ admits a representation of the form (fTJ) involving the following pair of Kraus operators, 



E = |0)(0| + VW|1)(1| 

E x = v^^|0)(l|, (7) 

By recourse to equation ([2]) the alluded Kraus representation of the single-qubit dynamics can be used to obtain the 
full dynamics characterizing the whole n-qubit system. 

III. TYPICAL ENTANGLEMENT DYNAMICS FOR n QUBITS 

A. Preliminaries 



We are going to study the average entanglement dynamics associated with random initial multiqubit pure states. 
These initial states arc of the form 



I*) = I7< n >|*„>, (8) 

where l^o) is a fixed n-qubit state, and U^ 1 ' is a random unitary transformation acting on the multiqubit system. The 
unitary transformations [/(") are generated according to the uniform distribution determined by the Haar measure 
(see [38[ and references therein). To evaluate the average properties of the entanglement dynamics of the n-qubit 
system we determine the entanglement evolutions corresponding to many initial states generated according to the 
procedure explained, and compute the concomitant mean entanglement evolution. 
We also investigate the entanglement dynamics of random mixed states of the form, 



/Wed = X|*)<*| + 2 „ I (9) 

with < x < 1. Here x is the purity of the multipartite state, I is the n x n identity operator, and |\&) is a random 
pure state generated following the same procedure as the one used to generate random initial pure states. 

To quantify the entanglement of a multipartite state one can deal with measures defined as the average of bipartite 
entanglement measures over all the possible bi-partitions of the system. This is usually known as global entanglement 
and a considerable amount of research has recently been devoted to its study (see [H, H(| and references therein) . 
The global entanglement E can be expressed as follow, 



1 [n/2] 

E = m Y. m 

L / J m=1 

E(m) = E ^(bipart.). (11) 

bipart bipart 

In equations (|10llll[) the possible bi-partitions of the n-qubit system are grouped in families corresponding to partitions 
into two sub-systems with given numbers of constituting qubits. ij( m ) denotes the average entanglement associated 
with the particular family comprising the K^^ art non-equivalent bi-partitions of the system into two subsystems 

having, respectively, m and n—m qubits. The sum appearing in equation (| 1 1 j) has, therefore, K^^ art terms. -E(bipart.) 
stands for the entanglement associated with one single bi-partition of the n-qubit system. 
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The above procedure for evaluating a global entanglement measure on the basis of bi-partitions can be implemented 
in a variety of ways, depending on which measure -E(bipart.) one uses to evaluate the entanglement of the individual 
bi-partitions. When dealing with mixed states the negativity constitutes the most important and most widely used, 
practical bipartite entanglement measure. The negativity iV(bipart.) corresponding a given bipartition of the system 
into two subsystems, one with m qubits and the other one with n — m qubits, with to < n — to, is given (up to 
a normalization multiplicative factor) by the sum of the negative eigenvalues at of the partial transpose matrix 
associated with this bipartition, 

tf(bipart.) = 2^EM- (12) 

i 

In this work we also consider an indicator of entanglement Em b proposed by Mintert and Buchleitner [41( which 
is given by 

E MB \p] = 2Tr[p 2 } - Tr[p\] - Tr[p% (13) 

with p\ and pi the reduced density matrices of the associated subsystems. The entanglement indicator Emb consti- 
tutes a lower bound of the squared concurrence, which for mixed states p is defined as [42j 

C[p]=miJ2piCm (14) 

i 

where {pi,^i} Q^i-P* = -0 are au the possible decompositions of p into mixtures of pure states and C is the 

concurrence of a pure state of a bipartite system, 

= y/2[l Tr{pl)] (15) 
whit p r the marginal density matrix corresponding to either subsystem. 

B. Pure states 

In this section we determine numerically the decay of the average entanglement of a multiqubit system when 
considering general, random initial pure states uniformly generated according to the Haar measure. All exhibited 
curves are obtained generating (and averaging over) 10 4 random pure states. See the Appendix A for technical 
details on this choice of the sample size. The main results obtained by us are summarized in figures 1-4. In figures 
la-Id we plot (for systems of three, four, five and six qubits respectively) the average entanglement (A) and the 
dispersion AA = yj (N 2 ) — (^V) 2 as a function of the dimensionless quantity jot. The dispersion AJV gives an 
estimation of the width of the entanglement distribution associated with initial random states. In all cases it is seen 
that AN is relatively small compared with (A), which means that this last quantity is representative of the typical 
entanglement dynamics exhibited by the multiqubit system. In figure 1 we also depict the curves corresponding to 
the entanglement decay of some important particular states, such as the n-qubit GHZ state (which is the maximally 
entangled state in the three qubit case) and the maximally entangled states for four, five, and six qubit systems 
(the corresponding curves in figures lc-ld are identified with the label "robust", for reasons explained below). For 
systems of four qubits we compute the entanglement evolution of the state HS, which was introduced by Higuchi and 
Sudbery in [43| (see also [4(1). This state has been conjectured to maximize the entanglement of four-qubit states 
[43| and, even if a proof of this conjecture is still lacking, various analytical and numerical evidences supporting it 
have been reported in the literature [H, 0, |43| . In the cases of five and six qubit-systems we plot the entanglement 
decay for the maximally entanglement states |$^f ) HU and |$^f ) [So[. The HS four-qubit state and the |$*^f ) six 
qubit states are specially relevant in connection with the process of decoherence because they also constitute robust 
states under decoherence. In a previous work [19( we have numerically determined these robust states for systems 
of four, five, and six qubits under different (Markovian) decoherence processes. By a "robust state" we mean an 
initial state evolving under decoherence to states that are, at any given time, more entangled than those generated 
by alternative initial states. The five qubit state I'&^f) does not coincide with the robust state reported in [li|, but 
in the case of the dynamics studied in the present work, our numerical results indicate that both states lead to the 
same (or very close) average entanglement evolution. 

A particularly interesting aspect of figures la-Id is the status of the state GHZ with respect to the typical 
behaviour of the entanglement dynamics of the multiqubit system. In the case of three-qubit systems the time 
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FIG. 1: Entanglement evolution of the 3, 4, 5 and 6 qubits random (pure) states under decoherence. All depicted quantities 
are dimensionless. 



dependent average entanglement is appreciably below the entanglement of states that evolved from the GHZ one, 
particularly at those times when the system's entanglement reaches its maximum values. The situation is different 
in the case of four-qubits, where we found that the entanglement dynamics associated with a GHZ initial state is 
quite representative of the typical, average dynamics associated with random initial states. On the other hand, for 
systems of five and six qubits the entanglement evolution associated with the GHZ state does not reflect the average 
behaviour, which is more closely represented by the entanglement dynamics exhibited by states respectively evolving 
from the previously mentioned l^j^jf ) and |3>g£^) states. 

Evidence has been previously reported suggesting that in some situations the robustness of entanglement may 



increase with the number of qubits in the system |8[. However, recent investigations 12[ indicate that, even when 
entanglement takes a long time to disappear, it rapidly adopts values too small to be of any practical relevance for 
technological applications. It is interesting to explore the dependence of the entanglement robustness on the number 
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FIG. 2: Entanglement evolution of the 3-qubit and 6-qubit representative states under decoherence: (a) GHZ states of three 
and six qubits; (b) random pure states of three and six qubits. All depicted quantities are dimensionless. 



of qubits within the context of the non-Markovian evolution studied in the present work. Some results pertaining 
to this issue are summarized in figures 2 and 3. In figure 2. a we compare the entanglement evolution yielded by 
initial GHZ states of systems of three and six qubits. It transpires from this figure that the finite time intervals 
of zero entanglement are longer in the case of six qubits than in the case of three qubits. This suggests that the 
lengths of these intervals tend to increase with the number of qubits, consistently with other results for particular 
states recently reported in [321 ] . This trend, however, is not observed when we evaluate the average entanglement 
evolution corresponding to random initial pure states, as can be appreciated in figure 2.b. Furthermore, the average 
entanglement evolutions corresponding to random initial pure states of systems consisting of 3, 4, 5, and 6 qubits are 
depicted in figure 3, showing an almost universal pattern. The average entanglement of random initial pure states 
thus seems to be less fragile than the entanglement corresponding to the GHZ state. It is interesting to compare 
these findings with those reported in a recent study by Aolita et al. [Hj]. These authors investigated the decay of 
the entanglement under various decoherence channels of the GHZ state. Aolita et al. determined an upper bound for 
the entanglement of states evolved (under the alluded decoherence channels) from initial GHZ states. They found, 
however, that the average entanglement of states evolved from random initial pure states does not comply with this 
bound, suggesting that the exponentially increasing fragility of the GHZ entanglement with the number of qubits of 
the system is not a generic feature of typical pure states. Our present results constitute further evidence also pointing 
towards the conclusion that the exponential (in terms of system's size) fragility of the GHZ state may be a peculiarity 
of this state not shared by generic ones. On the contrary, general pure states, on average, seem to exhibit more robust 
entanglement than the GHZ state. Another finding that is relevant in connection with this problem was reported 
by Hein, Diir and Briegel in [Tl| . where it was shown that the entanglement decay under decoherence corresponding 
to certain graph states states is very different from the one characterizing the GHZ states. The entanglement decay 
associated with graph states with a constant degree (not dependent on the number of qubits) does not become more 
acute as the size of the system increases. Results indicating that the entanglement fragility of the GHZ states is not 
a generic feature of multi-qubit states were also presented by Carvalho, Mintert and Buchleitner Q who showed that 
the rates of (exponential) decay of entanglement of the W sates are independent of the number of qubits of the system 
in the cases of coupling to a zero temperature bath and of dephasing. On the contrary, under these environmental 
influences the entanglement decay rate grows linearly with the number of qubits for the GHZ state. 

Let us denote by (N^) the average negativity corresponding to all the partitions of the system into two 
subsystems having, respectively, m and n — m qubits. In figured] we compare for 4 and 5 qubit systems the evolution 
of the entanglement of the most balanced bipartitions (N^) and the most unbalanced ones (iV^ 1 )). We can see that 
the behavior exhibited by the negativities corresponding to these special bipartitions is very similar to that of the 
global negativity (N) . 

The dynamical evolution of the average (Emb) is compared with that of the negativity in figure [5] for 4 and 5 
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FIG. 3: Average entanglement evolution under decoherence corresponding to random initial pure states of systems of 3,4,5, 
and 6 qubits. All depicted quantities are dimensionless. 



random qubit systems. In both cases the quantity Emb is able to detect entanglement only during the first period 
of entanglement decrease, but it does not detect the subsequent entanglement revivals. This is consistent with the 
fact that the quality of the entanglement indicator Emb tends to deteriorate as one considers states of increasing 



mixedness 4f| , and the multi-qubit systems under consideration here have a considerable degree of mixcdncss during 



the entanglement revival events. 



It is known that the entanglement exhibited by composite quantum systems tends to decrease as we consider 
states with increasing degree of mixedness. It is interesting to notice that the average entanglement evolution of the 
system under consideration in the present work does not follow this general tendency. The average time evolution of 
both the negativity and the degree of mixedness (as measured by the linear entropy Sl of the time dependent state) 
corresponding to random initial pure states of six qubits are depicted in figure 6. We can see that, except for the initial 
entanglement decrease, the moments of high entanglement tend to coincide with moments of high mixedness. This 
point was recently highlighted in [35], in connection with the evolution of systems of two-qubits, and is also consistent 
with previous results advanced in [32j | . Before discussing in more detail this aspect of the evolution of the multi-qubit 
systems discussed here it is worth mentioning that this behavior is atypical but should not be regarded as paradoxical. 
Even if, as a general tendency, entanglement tends to decrease when the degree of mixedeness increases, there is no 
strict correlation between these two quantities. It is indeed possible to find pairs of states such that the state with 
higher mixedness is also the more entangled one. In spite of this, it is still interesting that the system's evolution 
considered here consists (during certain time intervals) of states showing this somewhat atypical characteristic. 

To understand better this peculiarity let us first notice that the system's entanglement vanishes at those times 
when the parameter p t becomes equal to zero, pt = 0. At these times the mixedness of the multi-qubit system also 
vanishes, because its state is given by p s = (|0)(0|)® n , which is pure and separable. When starting from an initial 
entangled pure state pi = the evolution of the system's mixedness can be construed as arising from two 

opposite tendencies. First, due to the interaction with its surroundings, some entanglement is developed between 
the system and the environment and the system's mixedness tends to increase. On the other hand, as the evolution 
proceeds, the system eventually begins to approach the pure state p s , and its mixedness, after having achieved 
a maximum value, starts to decrease again. Consequently, the resulting evolution comprises an initial phase of 
monotonous increase in the degree of mixedness followed by a period of decreasing mixedness during which its value 
goes back to zero. On the other hand, during this whole time interval before pt achieves its first zero, the system's 
entanglement decreases monotonously, finally vanishing when p t = 0. Notice that during the second half of this 
period, during which the system is evolving towards p s , the system's entanglement and mixedness both decrease 
together and finally vanish when p t = 0. After this, when p t starts increasing again during the first entanglement 
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FIG. 4: Entanglement evolution of the random 4 (lines) and 5 (dots) qubit systems under decoherence for the most unbalanced 
(solid line, full dots) and most balanced (dotted line, empty dots) possible bipartitions of the systems. All depicted quantities 
are dimensionless. 



4-qubit random pure states 



5-qubit random pure states 




0.2 

o.o 

-0.2 
-0.4 

-0.6 





— <n> : 




<E > 




MB 





r * 



r * 



FIG. 5: The negativity and the Emb entanglement indicator as a function of 7oi for random 4 (a) and 5 (b) initially pure qubit 
systems. All depicted quantities are dimensionless. 



revival, the system's entanglement and degree of mixedness grow together, reaching their respective maximum values 
approximately at the same time and then decreasing together as the system evolves back to the state p s . 
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FIG. 6: Average time evolution of the global entanglement and the linear entropy corresponding to random initial 6-qubit 
pure states (a) and the mean global entanglement vs. the mean linear entropy associated with those states (b). All depicted 
quantities are dimensionless. 




FIG. 7: Time evolution of {N} for random initial n-qubit mixed states of the form Q as a function of 'yot and x, for n = 3 
(left) and n = 6 (right). All depicted quantities are dimensionless. 



C. Mixed Initial States 

Here we consider the entanglement decay of states of the form , consisting of a statistical mixture of a pure state 
| "J) and the totally mixed state. Our results are given in figure 7, where the average negativity (N) is plotted against 
7oi and the parameter x for systems of three and six qubits. We see that the general patterns associated with these 
two cases are similar except for the fact that the average negativity of systems having six qubits decreases more slowly 
with x than the average negativity corresponding to systems of three qubits. 

IV. CONCLUSIONS 

When considering the average behaviour of generic initial pure states of multiqubit systems for n — 3,4,5, and 
6 qubits one obtains in each case a very similar, almost universal pattern of entanglement evolution. This is an 
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interesting result, specially in connection with some previous studies suggesting that entanglement becomes more 
robust as the size of the system increases. It is clear that the system under study here does not follow that trend. 
Nor does the average entanglement become appreciably more fragile as the number of qubits increases, at least within 
the range of n- values considered by us. 

In the cases of systems having 5 and 6 qubits, the robust states studied in [l9| exhibit, under the dynamics 
considered in the present work, a behaviour that is quite close to the average behaviour corresponding to generic 
initial pure states. In this sense, the entanglement dynamics of the alluded robust states is representative of the 
typical dynamics of initial pure states. On the contrary, the average entanglement dynamics of 5 and 6 qubits departs 
from the one exhibited by the corresponding GHZ states, whose entanglement appears to be more fragile than the 
entanglement of typical, average, pure states. These findings are consistent with other recent results in the literature, 
also indicating that the entanglement fragility of GHZ states under various decoherence channels is not a generic 
feature of typical multi-qubit pure states. 

Systems of four qubits show some peculiarities. In the case of these systems the average entanglement evolution of 
random initial pure states is quite similar to the entanglement evolution corresponding to an initial GHZ state. 

The observed universal features exhibited by the average entanglement dynamics associated with random initial 
pure states of n-qubits is plausibly related to the fact that each qubits is interacting with its "own" independent heat 
bath. This fact implies that the typical behaviours of systems with different (total) numbers n\ and rii of qubits, 
even though they are not identical, share some important features: the evolution rules for the local states of blocks 
of L qubits (with L < min(rii, 712)) are the same in both cases. For instance, the evolution equation for single qubits 
is always given by equation ([1]). 

The fact that the GHZ states deviate from the typical behaviour, having entanglement features that are more 
fragile than those of typical initial states, is consistent with other known manifestations of the "fragility" of the 
quantum correlations exhibited by these states. For instance, it is well-known that if we trace over one of the qubits 
of a system initially described by a GHZ state, the remaining n — 1 qubits are left in a mixed non-entangled state. 
Also, if we measure the state of one of the qubits, the remaining qubits are left in a non-entangled state. These 
properties are in contrast with the ones exhibited by other multi-qubit entangled states, such as the W states, 
where some entanglement survives after measuring the state of one of the parts of the system. As a final remark 
concerning GHZ states, it is worth mentioning that for systems having five or more qubits the difference between the 
entanglement of the maximally entangled states (according to the partitions-based measures considered here) and 
the GHZ states is much larger than in the case of systems of four qubits. This may be related to the fact that the 
behaviour of GHZ states deviate more from the average behaviour in the case of systems with five or more qubits, 
than in the case of systems with four qubits. 

The main results presented in this work have been obtained using entanglement measures based upon the 
negativities associated with the different possible bi-partitions of the multi-qubit system. We focused mostly 
on the behavior exhibited by the global mean value corresponding to all those bi-partitions. However, the en- 
tanglement dynamics corresponding to particular bi-partitions was also considered, and the results obtained are 
similar to those corresponding to the global average. The negativity, even though it provides one particular way 
of quantitatively characterizing the amount of entanglement in composite quantum systems, it constitutes one of 
the most widely used practical tools for studying the entanglement of mixed states. It would be interesting to 
explore other ways of characterizing the evolution of the entanglement (or other manifestations of quantum me- 
chanical correlations) of the system considered here. Any further developments along these lines will be very wellcome. 
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FIG. 8: {N) as a function of the sample size for 4 qubit random states. Inset: (N) as a function of the sample size for 3 (circle), 
5 (triangle), and 6 (square) qubit random states. All depicted quantities are dimensionless. 



Appendix A: Details of the Numerical Calculations 



In this appendix we provide some technical details on the numerical Monte Carlo procedure employed to compute 
the average properties corresponding to general initial pure states of the multi-qubit system. As was stated in the 
main text of this work, we compute these averages by generating a random sample of 10 4 initial pure sates. These 
states are generated according to the Haar measure (see [38| and references therein). The Haar measure determines 
a uniform distribution in the space of 2 n x 2" unitary matrices which, in turn, induces a natural uniform distribution 
on the Hilbert 2"-dimensional space corresponding to the n-qubit system. To generate uniformly distributed random 
initial states we generate random 2" x 2™ unitary matrices uniformly distributed according to the Haar 

measure, and then compute = U^ n '\^o), where |^o) is a Axed n-qubit state. 

To ensure that we have an appropriate sample, we consider increasing values of its size until the computed averages 
do not change any more (within a certain range of numerical error) with further increases of the sample size. The 
behavior of the numerically determined average negativity (N) as a function of the sample size, corresponding to the 
case of four qubits and 7oi = 40, is depicted as an illustration in figure 8. The behavior of (N) against the sample 
size corresponding to systems of three, five and six qubits is plotted in the inset of figure 8. It transpires from figure 
8 that the changes on the numerical values of (N) corresponding to an increase of the sample size from 10 4 to 10 5 are 
already negligible within the scale of the figures exhibited in the main text of this paper. 
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